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I.  Introduction 

Since  the  discovery  of  a large  dropout  (ionospheric  hole)  in  the 
ionospheric  total  electron  content  during  the  launch  of  NASA's  Sky lab 
(Mendillo  et  al.,  1975  a,b)  much  investigation  has  gone  into  the 
controlled  artificial  production  of  such  holes  or  plasma  density 
depletions  in  the  F region  ionosphere  (Bernhardt  et  al.,  1975; 
Bernhardt,  1976;  Bernhardt  and  da  Rosa,  1977;  Mendillo  and  Forbes, 
1978;  Pongratz  et  al.,  1978;  see  also  Trans.  Am.  Geophys.  Un. , 58, 
454-455,  1977,  Spring  Meeting  Special  Session  on  "Artificially 
Created  Holes  in  the  Ionosphere").  Typically  investigators  have 
focused  on  the  chemical  release  of  ^0,  ^ or  high  explosives.  These 
contain  molecular  species  which  react  with  ionospheric  0+  in  such  a 
manner  that  the  0+  loss  rate  due  to  these  reactions  is  several  orders 
of  magnitude  faster  than  natural  ionospheric  loss  rates  for  0+ 
(Mendillo  et  al.,  1975  a,b;  Bernhardt,  1976).  Recently,  Pongratz 
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et  al.  [1978]  conducted  two  high  explosive  (88kgra  mixture  of  nitro-  ,,s^~ 


methane  and  ammonium  nitrate)  F region  chemical  releases  (260  and  . i 

a m#  □ 

280  km  altitude),  codenamed  Lagopedo,  over  the  Hawaiian  Islands  Q 

(midlatitude)  and  formed  a hole  ("v  100  km  in  diameter)  via  the 


predicted  chemical  processes. 
Note:  Manuscript  submitted  March  22,  1978. 


Within  the  context  of  artificial  hole  formation  much  of  the 
theoretical  work  to  date  (Mendillo  et  al.t  1975  b;  Bernhardt,  1976; 
Pongratz  et  al.,  1978)  has  centered  on  the  complicated  chemical 
reactions  and  molecular  diffusion  equations  for  producing  these  holes 
(e.g.,  electrodynamic  forces,  such  as  polarization  electric  fields, 
on  the  ions  or  electrons  are  not  considered).  lu  the  past  few  years, 
experimental  (Kelly  et  al.,  1976;  McClure  et  al.,1977)  and  numerical 
simulation  studies  (Scannapieco  and  Ossakow,  1976)  have  observed 
(under  appropriate  conditions)  the  formation  of  natural,  rising 
plasma  density  depletions  (holes  or  bubbles)  in  the  equatorial 
Spread  F ionosphere.  In  the  present  paper  we  wish  to  perform  a 
marriage  between  the  two  studies  and  examine  the  behavior  of  a large 
plasma  density  depletion  (hole)  artificially  injected  into  the  bottom- 
side  nighttime  equatorial  F region  ionosphere. 

II . Theoretical  Model  and  Numerical  Simulation  Studies 

In  the  present  simplified  model  we  will  take  the  approach  that 
the  large  hole  has  already  been  created  and  thus  represents  a new 
Initial  condition  for  our  previous  nonlinear  numerical  simulation 
equatorial  Spread  F studies  [Scannapieco  and  Ossakow,  1976;  Ossakow 
et  al.,  1978  a,bj.  In  this  way  we  do  not  follow  the  chemical  kinetics 
of  the  actual  hole  making  process,  rather  we  follow  the  evolution  of 
the  hole  after  it  has  been  produced.  Pongratz  et  al.  [1978)  noted 
that  the  entire  hole  in  Lagopedo  was  produced  in  less  than  5 minutes 
and  Sjolander  et  al.  [1978]  noted,  using  rocket  in  situ  ion  mass 
spectrometer  measurements  during  Lagopedo,  that  the  dominant  ion  in 
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the  artificially  created  hole  was  0+.  Thus,  we  are  able  to  utilize 
the  two  dimensional  Cperpendicular  to  the  magnetic  field,  B)  numerical 
simulation  code  developed  for  natural  equatorial  Spread  F studies  in 
the  collisional  Rayleigh-Taylor  regime,  which  is  represented  by  the 
following  equations  (see  Scannapieco  and  Ossakow,  197b;  Ossakow  et  al., 
1978  a.b) 

f?  " | ^l  * *>  *7n  " -'Vn'‘no)  <]> 

V-Cv^nV^p  - ® (g  x z)-Vn  (2) 

where  n is  the  electron  density  (no  refers  to  the  background  equilib- 
rium density;  n = nQ  + n^),  is  the  polorization  (induced)  potential, 
vR  is  the  usual  (used  in  natural  Spread  F studies)  ionospheric 
recombination  rate  (see  Fig.  1),  vln  is  the  usual  ion-neutral 
collision  frequency,  g is  gravity  (pointing  in  the  -y  direction), 
z - B/|B|  and  V refers  to  the  x and  y directions.  In  our  simulations 

realistic  altitude  profiles  of  v and  v (see  Fig.  1)  and  n (see 

in  K o 

Fig.  2)  have  been  used. 

The  numerical  simulation  plane  coincided  with  the  equatorial 
plane  spanning  an  altitude  (y)  range  from  172  to  452  km  in  one  case 
and  252  to  532  km  in  the  other  case  (both  with  Ay  - 2 km)  with  an 
east-west  (x)  extent  of  200  km  (Ax  ■ 5 km).  The  system  of  eqns.  (1) 
and  (2)  was  initialized  as  follows 


• V -^^Mz:-JZi2Sf:.rT32:  T^g^T*****  y* 


n/n  “ 
o 


0.03  , r-rQ  < 13  km 

expIa(r-rQ-15)^]  - 0.97  , 15  km  < r-rQ  < 35  km 

1.0,  r-r  > 35  km  (3) 

o 
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where  a - ln(1.97)/(20)  . This  provided  a constant  97%  depletion 
over  a central  30  km  of  the  hole  and  the  depletion  decreased  to 
ambient  density  over  a 20  km  extent  on  each  side  (total  extent  of 
bubble  was  70  km).  The  thickness  of  the  hole  at  half  depletion  was 
"v  63  km  (see  for  example  Pongratz  et  al.,  1978).  The  results  of  a 
simulation  with  rQ  corresponding  to  xq  * 0,  yQ  “ 300  km,  i.e.,  release 
altitude  = 300  km,  are  shown  in  Figs.  2-5,  which  depict  contour  plots 
of  n/nQ  (or  equivalently  n^/no)  at  t * 0,  500,  1000,  and  2000  sec. 

The  t *•  0 frame  shows  the  initial  state  after  the  hole  has  been 
formed  by  chemical  processes.  At  t * 500  sec  the  hole  has  begun  to 
steepen  on  its  topside  as  it  rises  toward  the  F peak  (altitude  * 

354  km).  This  corresponds  to  a rise  velocity  *v»  28  m/s ec.  We  note 
that  the  initial  hole  is  no  longer  circular  and  plasma  density 
enhancements  are  forming  to  the  sides  of  the  bubble  (hole) . As  tne 
initial  depletion  rises  on  the  bottomside  it  displaces  regions  of 
higher  density  which  fall  to  lower  altitudes  creating  enhancements. 

At  t ■ 1000  sec  the  top  of  the  bubble  is  at  the  F peak  and  the  enhance- 
ments are  getting  larger.  This  should  be  contrasted  with  the  natural 
Spread  F developing  bubble,  in  this  ionospheric  profile,  which  took 
"v  8000  sec  to  reach  the  peak  (see  Scann&pieco  and  Ossakow,  1976; 

Ossakow  et  al.,  1978b).  At  t ■ 2000  sec  the  hole  is  well  past  the 
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peak  and  has  bifurcated  on  ita  topaide.  This  is  similar  to  hariuiu 
cloud  striation  bifurcation  phenomena  (aee  Zabuaky  at  al..  197 3; 
Ossakow  et  al.,  1977)  and  would  lead  to  smaller  scale  struaturea 
if  the  present  simulation  had  more  resolution.  The  hole  is  now 
accompanied  by  large  enhancements  and  the  enhancements  in  turn  are 
forcing  depletions  to  lower  altitudes.  This  picture  is  tar  different 
from  the  initial  circular  hole  at  t ■ 0.  This  figure  shows  a wide- 
spread rapid  equatorial  Spread  F condition. 

Figures  6-9  exhibit  the  results  for  an  artificial  depletion 
released  at  380  km  altitude  (x^  “ O,  - 180  km)  in  an  ionosphere 
with  the  altitude  of  the  F peak  equal  to  434  km.  in  this  case  the 
time  evolution  of  the  hole  is  much  more  rapid.  For  example, between 
t • 0 and  50  sec  the  hole  is  rising  with  a velocity  v 88  m/sec.  The 
initial  rise  velocity,  here  and  in  the  previous  case,  is  consistent 
with  a bubble  rise  velocity  given  by  (Oasakow  and  Chaturvedl , 1978) 

Vg  " (g/v^)  f (n^/nQ),  where  f (itj/n^)  Is  an  Increasing  function  oi 
the  percentage  depletion,  n^/n^,  and  dependent  on  bubble  shape.  In 
particular  f corresponds  here  to  a circular- 1 ike  shape.  The  case 
depicted  in  Figs.  6-9  is  also  to  be  contrasted  to  the  natural  Spread  F 
case  in  this  geometry  where  a bubble  took  ' 1000  sec  to  reach  the 
F peak. 

III.  Concluding  Remarks 

We  have  presented  numerical  simulation  results,  in  the 
coll isional  Rayleigh-Tay lor  regime,  of  artificial  depletions  injected 
into  the  bottomside  nighttime  equatorial  F region.  It  has  been  shown 
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that  such  depletions,  while  they  steepen  and  rise,  will  cause  plasma 
density  enhancements  to  be  formed.  The  entire  final  picture  is 
similar  to  natural  equatorial  Spread  F numerical  simulation  phenomena 
except  that  by  artificial  injection  it  occurs  on  a much  faster  time 
scale. 
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Fig.  2 Contour  plot*  of  constant  n(/n0  nt  1 **  0 sec  (initial  condition).  The 
circular  contours  represent  the  Artificial  depletion  (-)  Midi  thnt  the  outermost  con- 
tour m n 16%  depletion  and  the  innermost  is  n 96%  depletion  (inaidc  this  contour 
la  n constant  97%  depletion).  The  large  dashed  curve  represents  a plot  of  the  am- 
bient electron  number  density  (values  on  upper  horizontal  axis),  n„,  as  a function 
of  altitude.  The  vertical  (y)  axis  represents  altitude,  the  lower  horizontal  (x)  axis 
east-west  range,  and  the  ambient  magnetic  field  is  out  of  the  figure  (z). 
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Ki)j.  3 Sumo  ns  Ki*|.  (1).  except  t - 500  «v,  The  smaller  dashed  contours 
are  contours  of  constant  n,/n„  for  plasm*  density  enhancements  (•*■).  The 
enhancement  contour  repn*senU  a 111%  enhancement. 
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